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ABSTRACT 

We present dynamical modeling of the broad line region (BLR) in the Seyfert 1 galaxy Mrk 50 
using reverberation mapping data taken as part of the Lick AGN Monitoring Project (LAMP) 2011. 
We model the reverberation mapping data directly, constraining the geometry and kinematics of the 
BLR, as well as deriving a black hole mass estimate that does not depend on a normalizing factor or 
virial coefficient. We find that the geometry of the BLR in Mrk 50 is a nearly face-on thick disk, with 
a mean radius of 9.6^'g light days, a width of the BLR of 6.9^'^ light days, and a disk opening angle 
of 25 ± 10 degrees above the plane. We also constrain the inclination angle to be degrees, close 



to face-on. Finally, the black hole mass of Mrk 50 is inferred to be log 10 (MBH/^o) = 7.57 



+0.44 
-0.27- 



By 



comparison to the virial black hole mass estimate from traditional reverberation mapping analysis, 
we find the normalizing constant (virial coefficient) to be log 10 / = 0.78^° 27) consistent with the 
commonly adopted mean value of 0.74 based on aligning the Mbh-c* relation for AGN and quiescent 
galaxies. While our dynamical model includes the possibility of a net inflow or outflow in the BLR, 
we cannot distinguish between these two scenarios. 

Subject headings: galaxies: active — galaxies: individual (Mrk 50) — galaxies: nuclei 
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The standard model of active galactic nuclei (AGNs; 
lAntonuccl 119931 : lUrrv fc Padovanil Il995| ) explains their 
broad emission lines as being produced in a broad emis- 
sion line region (BLR) situated on the order of light days 
from the black hole dWandel et al.l Il999t iKaspi et al.l 
120001: iBentz et al.l l2006f h" The distance of the BLR 
from the black hole can be measured using reverbera- 
tion mapping, in which the average delay time is mea- 
sured between a timeseries of the variable AGN contin- 
uum luminosi ty and a timeseries of the variable broad 
line emission dBlandford fc McKedfl98l !Petersodll99l: 
iPeterson et al.ll2004D . Standard reverberation mapping 
analysis also provides estimates of the black hole mass, 
Mbh, in AGNs from a normalized virial product. The 
virial product, M v i r = / v 2 ct /G, is derived from the 
average light travel time lag of the BLR, r, and the 
typical velocity of the BLR gas, v, measured from the 
width of the broad lines. The small sample of ~ 50 
reverberation mapped AGNs is then used to deter- 
mine single-epoch A/bh estimates for much larger sam- 
ples of AGNs using the BLR-size to luminosity rela- 
tion (IVestergaard fc Peterson! 120061: iMcGifl et ail 120081: 
IVestergaardlboilh . 

However, there are certain limitations to the stan- 
dard reverberation mapping techniques. The object- 
to-object scatter of the normalizat ion factor / is be- 



lieved to be of order ~ 0.4 dex (lOnken et al 



Collin et al.ll20M I Woo et alJl2010t [Greene et al.1 



[2004 : 



2010a: 



Graham et alJl20*TT ) based on assumin g the same M] 



BH" 



relation (e.g. iBennert et all 1201 lh as for quiescent 
galaxies. It would be desirable to avoid this assump- 
tion and estimate -Mbh from reverberation mapping data 
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alone. The details of the BLR geometry and dynamics 
are also poorly constrained by standard reverberation 
mapping analysis. Measuring the time lag as a func- 
tion of line-of-sight velocity has shown that while some 
BLRs are consistent with virial mo t ion in a Keplerian 
potential (iPeterson fc Wandell 119991 : iBentz et all 120091 : 
iDennev et al.l I2010D some show suggestions o f inflow- 
ing gas (|Bentz et al.|[2010t IDennev et al.ll2010t) . In ad- 
dition to the mean radius of the BLR as obtained in the 
standard analysis, we would like to constrain the overall 
geometry of the BLR in more detail. 

Recent improvements in reverberation mapping data 
and analysis are starting to provide better constraints 
on the geometry and dynamics of the BLR. Velocity- 
resolved transfer functions (VRTFs) have been mea- 
sured using high-quality reverberation mapping data 
from the Lick AGN Monitoring Project in 2008 (LAMP 
2008. IWalsh et al.ll2009HBentz et alj|2009l) and from the 
2007 MDM Observatory reverberation mapping cam- 
paign (jDennev et al.1 l2010T l , showing signatures consis- 
tent with inflow, outflow, and virialized motion for 
different AGNs. However, a clear interpretation of 
VRTFs requires additional modeling steps, since they 
are functions of time lag instead of position within the 
BLR. The traditional reverberation mapping a nalysi s 
has also been recently improved by IZu et al.l (|2011l ). 
who model the AGN continuum and line light curve us- 
ing an implementation equivalent to Gaussian Processes 
(iKellv et al1l2q09UKozlowski et alJl2010b iMacLeod et al.1 
120101 IZu et al.l I2012T ). Members of our team have de- 
veloped a method for determining the geometry and dy- 
namics of the BLR by directly modeling reverberation 
mapping data (jPancoast et alj|2011t iBrewer et aIll201lL 
hereafter Pll and Bll respectively), estimating the un- 
certainties in the framework of Baycsian statistics. Our 
modeling method constrains -Mbh without requiring a 
normalization constant /. We also constrain the geom- 
etry of the BLR, its orientation with respect to the line 
of sight, and the possibility of net inflowing or outflow- 
ing gas in the BLR. We have previously demonstrated 
our method on LAMP 2008 data for Arp 151 and esti- 
mated A/bh with smaller uncertainties than traditional 
reverberation mapping analysis (Bll). 

What is now needed to make further progress is large 
samples of high quality velocity resolved reverberation 
mapping data. For this purpose we carried out an 11- 
week spectroscopic observing campaign at Lick Obser- 
vatory, the Lick AGN Monitoring Project 2011. The 
project focused on nearby AGNs with bright H/3 lines, 
which are good candidates for dynamical modeling. Here 
we present the first results of dynamical modeling for the 
project, focusing on one of the most variable objects in 
the sample, Mrk 50. The average time lag and virial 
Mbh estimates from tra ditional reverberation mapping 
analysis are presented by iBarth et al.1 (|2011bl ) . Here we 
present an alternative analysis based on our direct mod- 
eling technique. The H/3 and V^-band continuum light 
curve data are briefly described in Section 2, the dynam- 
ical model for the BLR is described in Section 3, and our 
results and conclusions are given in Section 4. 

2. DATA 

We observed Mrk 50 in Spring 2011. The data, shown 
in Figurc[T] and the top left panel of Figurc[2j include a 
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Fig. 1. — The integrated H/3 broad line and V-band con- 
tinuum light curves. The H/3 light curve has flux units of 
10 — 15 ergem - 2 s . The V-band light curve is in arbitrary flux 
units. 



light curve of V^-band continuum flux and a time series 
of the broad H/3 line spectral profile. More observational 
details, as well as details about the measur ement of V- 
band a nd H/3 light curves, are described bv lBarth et al.l 
(|2011bf ). We model all 156 epochs of the F-band light 
curve and 55 of the 59 epochs of the H/3 line profile, 
ignoring those epochs with low S/N or other problems. 
The median S/N for the H/3 line profile throughout the 
campaign is 75 per pixel. 

AGN and stellar continuum lines can significantly alter 
the measured broad line widths in AGNs, affecting sin- 
gle ep och .Mbh estimates (jDennev et al.| [2009; Par k et all 
|2011|) . In order to reduce contamination from other lines 
when modeling the H/3 line profiles, the Mrk 50 spec- 
tra have been fit with AGN and stellar continuum com- 
ponents and the He II A4686 line just blueward of H/3, 
and then these co mponents were subtr acted to yield the 
"pure" H/3 profile (|Barth et al.ll201lbn . 

3. THE DYNAMICAL MODEL OF THE BLR 

We now give a brief description of our method for di- 
rectly modeling reverberation mapping data. The mo- 
tivation for our approach is developed in Pll and fur- 
ther implementation details are described in Bll. We 
model the BLR as a large number of point-like clouds, 
each with a given position and velocity. Several parame- 
ters describe the overall spatial distribution of the clouds 
and the prescription for assigning velocities to the clouds, 
given their positions. Our goal is to estimate these pa- 
rameters. 

The continuum emission from the central ionizing 
source is absorbed by these clouds and re-emitted as 
broad line emission, allowing us to predict the line flux 
and shape as a function of time, i.e. to produce mock 
data sets of the form shown in Figure [21 

The full set of model parameters includes the geometry 
and dynamics parameters for the BLR clouds corrected 
to the rest frame of Mrk 50, as well as a continuous ver- 
sion of the continuum light curve, since the continuum 
light curve must be evaluated at arbitrary times in order 
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Fig. 2. — Top: H/3 spectra in velocity units for each epoch in the 
light curve for data, left panel, and model, right panel. Dark red 
corresponds to the highest levels of flux and dark blue corresponds 
to the lowest levels, where the same color scale is used for the data 
and model. Middle: integrated H/3 flux for each epoch in the light 
curve for the data, blue solid line with errorbars, and model, red 
dashed line. As an illustration of the range of solutions, we show 
light curves for five acceptable models as dotted gray lines. For the 
correct time separation between light curve epochs, see Figure [T] 
The model is able to reproduce the major features of the data. 
Bottom: two examples of H/3 spectra fit by the model, with data 
shown by blue and green errorbars and model fits shown by red 
lines. 
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Fig. 3. — Geometry of the BLR for three models, with the x, y, 
and z axis scales in light days and the observer's line of sight along 
the x-axis. The top panel BLR distribution is a close to face-on 
torus of clouds, the middle BLR distribution is a close to face-on 
disk of clouds similar to the geometry inferred for Mrk 50, and the 
bottom BLR distribution is a dense sphere of clouds. 



to compute mock data for comparison with the actual 
data. The observed continuum light curve is interpo- 
lated using Gaussian Processes to create a continuous 
light curve and to account for the uncertainty in the in- 
terpolation. Gaussian Processes have been found to be 
a good model for larger samples of AGN light curves 
dKellv et alJl2009HKoziowski et aIJl20Tot iMacLeod et al.l 
l2010HZu et alJl201lLl20l3T 

The model for the BLR geometry is simple yet flex- 
ible, allowing for disk-like or spherical geometries with 
asymmetric illumination of the gas. Examples of possi- 
ble BLR geometries are shown in Figure [3] The model 
for the spatial distribution of the BLR gas is first gener- 
ated from an axisymmetric 2-D configuration in the x-y 
plane, with a parameterized radial profile. The radius r 
of a cloud from the origin is generated as follows. First, 
a variable g is drawn from a Gamma distribution with 
shape parameter a and scale parameter 1: 

S~r(a,l) (1) 

Then, the radius r of the cloud is computed by applying 



the following linear transformation to g: 



r = F/J.+ 



m(i-f) 



(2) 



The parameters {/x, F, a} control the radial profile of the 
BLR. fj, is the overall mean radius of the BLR (this can 
be verified by taking the expectation value of r in Equa- 
tion^. The parameter F E [0,1] allows for the pos- 
sible existence of a hard lower limit Ffi on radius, be- 
cause there may be some radius interior to which the 
BLR gas would all be ionized and thus u nable to respond 
to ch anges in the continuum emission (|Korista fc GoacH 
12004 ). a controls the shape of the gamma distribution: a 
value of a close to 1 imposes an exponential distribution 
(allowing for disk or ball configurations) , whereas large 
values of a create a narrow normal distribution (allowing 
for shell or ring configurations). In the implementation, 
and in the description of the same model in Bll, we pa- 
rameterise the shape by /3 = 1/ 'y/a instead of a because 
P has a simple interpretation as the standard deviation 
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Fig. 4. — Inferred posterior PDFs for model parameters, including Mbhj inclination angle (0 degrees is face-on), and opening angle of 
the BLR disk. Joint posterior PDFs are also shown to illustrate the major degeneracies. 



of g in units of its mean. The radial width of the BLR 
can be defined as the standard deviation of r: 



G r = Hp(l-F). 



(3) 



In order to assign velocities to the BLR gas clouds, 
the model uses probabilistic perturbations about circular 
orbits. The solution for the radial velocity of a BLR cloud 
given its position r, energy E, and angular momentum 
L is: 




GM 



^BH 



£2 



(4) 



If we wished to impose circular orbits, the values for E 
and L would be fully determined by the radius r of the 



cloud: 




:±n 2 E 



(5) 
(O 



To obtain elliptical orbits, we generate values for E and 
L probabilistically, given r. The probability distributions 
for energy and angular momentum are parameterized by 
the parameter A and are given by: 



E = 



1 



1 + exp(-x) 



p(L) oc exp ( — 



E r 



L\ 



(7) 
(8) 
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Fig. 5. — Inferred posterior PDFs for model parameters, includ- 
ing the mean radius of the BLR, radial width of the BLR, and the 
inflow fraction of BLR gas. 
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Fig. 6. — Examples of acceptable transfer functions for Mrk 
50. The top two and middle left panels show examples of VRTFs 
drawn from the model parameter posterior PDFs, illustrating the 
range in inferred transfer function shape. In the color-code of the 
VRTFs, red corresponds to the highest levels of response and dark 
blue corresponds to the lowest levels. The middle right panel shows 
the mean lag fo r each of the VRTFs . The mean lag in seven ve- 
locity bins from[Barth et al. ( 201 lbl) are shown by red errorbars, 
which were measured by cross-correlation analysis. We calcu late 
the mean lag in the seven velocity bins of lBarth et~aLl j2011bl) for 
~ 200 VRTFs made using model parameters drawn randomly from 
their posterior PDFs, shown in light blue. The bottom panel shows 
the velocity-integrated transfer functions for the VRTFs shown in 
the first three panels. 



where \ ~ W(0,A 2 ) and \L\ < \L cil - c \. For A -> we 
recover circular orbits and increasing A creates more el- 
liptical orbits. Since there are two solutions for the sign 
of v r , the model also includes a parameter for the frac- 
tion of outflowing vs inflowing gas. The inflowing and 
outflowing gas is bound to the gravitational potential of 
the black hole, but an inequality in the fraction of inflow- 
ing and outflowing gas has the desired effect of modeling 
asymmetries in the H/3 spectral line profile as observed in 
Arp 151 (Bll) when an asymmetric illumination model 
is included. 

Once a 2-D configuration of clouds in the x-y plane has 
been generated, and velocities have been assigned to the 
clouds, rotations are applied to "puff up" the 2-D con- 
figuration into a 3-D configuration. We first rotate the 
cloud positions about the y axis by a small random an- 
gle; the typical size of these angles determines the open- 
ing angle of the cloud distribution. The opening angle 
is defined as the angle above the midplane of the disk 
or sphere. We then rotate around the z axis by random 
angles to restore the axisymmetry of the model. Finally, 
we rotate again about the y axis, by the inclination angle 
(common to all of the clouds) to model the inclination of 
the system with respect to the line of sight. The inclina- 



tion angle is defined so that zero degrees corresponds to 
a face-on configuration and 90 degrees corresponds to an 
edge-on configuration. 

In order to produce asymmetric broad line profiles, we 
include a simple prescription for asymmetric illumination 
of the BLR clouds. We assign a weight w to each cloud, 
given by w = 0.5 + k cos 4>, where 4> is the azimuthal po- 
sition of the cloud in spherical polar coordinates. The 
parameter k ranges from 0.5, corresponding to illumi- 
nating the near side of the BLR, to —0.5, corresponding 
to illuminating the far side of the BLR. Physically, the 
near side of the BLR could be preferentially illuminated 
if the far side of the BLR were obscured by gas, and the 
far side of the BLR could be preferentially illuminated if 
the BLR clouds only reradiate the continuum emission 
towards the central ionizing source due to self-shielding 
within the cloud. Inflowing gas with the near side of the 
BLR illuminated can, in principle, be distinguished from 
outflowing gas with the far side of the BLR illuminated 
by the VRTF, since the lags for these two cases are dif- 
ferent. 

In addition, we allow for a scaling factor to describe the 



6 



Pancoast et al. 



percentage variability of the emission line compared to 
that of the continuum. While for Arp 151 the variability 
of the continuum was approximately equal to that of the 
H/3 flux, in the case of Mrk 50 we find that the continuum 
variability is less than that of the line. This is consistent 
with the amplitude of variability of the ionizing contin- 
uum res ponsible for H/3 being larger than that of the 
V-band (jMeusinger et al.|[201lL and references therein). 

Once the dynamical model has been defined, we are 
able to compute simulated data that are then blurred 
with a Gaussian kernel to model the instrumental reso- 
lution. The simulated data are then compared with the 
actual data. For the likelihood function, we use the stan- 
dard Gaussian assumption: 



H. 



4. RESULTS AND CONCLUSIONS 



P(data|modcl) cx exp 



1 



X 2 (model, data) 



(9) 



With the likelihood function defined, the modeling 
problem is reduced to computing the inferences on 
all of the model parameters. The likelihood func- 
tion, P(data|model), is combined with the prior dis- 
tribution for the parameters using Bayes' Theorem: 
P(model|data) cx P(model) x P(data|model). The 
posterior probability distribution for the parameters is 
sampled using the Diffusive Neste d Sampling algorithm 
(|Brewer. Partav. fc Csanvlll2011b[) . Nested Sampling al- 
gorithms initially sample the prior distribution, and sub- 
sequently create and sample more constrained distribu- 
tions, climbing higher in likelihood. In the specific case 
of Diffusive Nested Sampling, uphill and downhill moves 
are allowed, allowing the exploring particles to return to 
the prior, take large steps, and then climb the likelihood 
function again. We assigned uniform priors to most pa- 
rameters except for the mean radius and Mbh, which 
have log uniform priors to describe initial uncertainty 
about the order of magnitude of the parameter. 

By computational necessity, our model is relatively 
simple. While it is still rather flexible and can repro- 
duce the large scale features of the reverberation map- 
ping data, it is unable to model every detail of the H/3 
light curve. The large scale features of the variability in 
the H/3 light curve are well-modeled, for example, but the 
small epoch-to-epoch fluctuations in the light curve are 
not (see Figure [2J. In addition, the errorbars reported on 
the data are very small, and our model is not able to fit 
the data set to within these small error bars (i.e. we can- 
not achieve reduced \ 2 ~ !)■ If we did not take this into 
account our uncertainties would be unrealistically small. 
This issue is a generic feature of the fitting of simply pa- 
rameterized models to informative data sets, and will be 
discussed in depth in a forthcoming contribution (Brewer 
et al, in preparation). In order to account for this effect 
and to obtain realistic and conservative uncertainties, we 
explore the effect of inflating the errorbars on the spec- 
trum data by a factor H, or equivalently, choosing to 
form the posterior distribution from different chunks of 
the Nested Sampling run (i.e. different ranges of allowed 
likelihood values) . For each value of H tested, we inspect 
the posterior distribution over simulated data (top right 
panel in Figure [5]) to ensure that the major features of 
the data are reproduced. We find that, as long as H is 
low enough that the models fit the major features in the 
data, the resulting posterior distributions on the param- 
eters are insensitive to the exact choice of the value for 



Our inferred geometry and dynamics parameters of the 
BLR in Mrk 50 are shown in Figures [4] and [5j The 
shape of the BLR gas radial profile is constrained to 
be closer to exponential (a < 1), with a mean radius 
of fi = 9.6^J'g light days and a width of a r = 6.9^^ i 
light days (the uncertainties quoted are symmetric 68% 
confidence limits). Even though the mean radius is not 
simply c times the mean lag in the general asymmetric 
case, we expect our mean radius to roughly correspond 
to the lag measuremen ts using cross-correlation analysis 
bv lBarth etall (|2011bf ). which arc r pcak = 9.751°;^ and 

Teen = 10.64^0.93 hght days. Our mean radius agrees 
more closely with r pca k, although r con is more commonly 
used for black hole mass estimation. We infer the inner 
radius of the BLR distribution to be Ffj, = light 
days. The opening angle of the BLR disk, defined be- 
tween 0-90 degrees, is 25 ± 10 degrees, closer to a thin 
disk than to a sphere. The inclination angle of the thick 
BLR disk with respect to the line of sight is constrained 
to be 9^5 degrees, closer to face-on, consistent with the 
standard model of bro ad line AGNs ()Antonuccil 119931 : 
lUrrv fc Padovanilll995[) . 

The dynamical modeling results constrain Mrk 50 to 
have 39% probability of net inflowing gas and 61% proba- 
bility of net outflowing gas, with equal amounts of inflow- 
ing and outflowing gas ruled out (inflow fraction = 0.5), 
as shown in Figure [5] This result suggests only a slight 
preference for outflow while the need for either outflow 
or inflow is quite robust, suggesting that a more phys- 
ical model for inflow and outflow is needed in order to 
distinguish between them for the case of Mrk 50. Equal 
amounts of inflowing and outflowing gas are ruled out 
because net inflowing or outflowing gas, along with the 
illumination model, creates the asymmetry in the H/3 line 
profile observed in the data. 

In addition to constraining the geometry of the BLR, 
our dynamical model also places an independent estimate 



+0.44 
-0.27- 



on M B h, inferred to be log 10 (M B H/M Q ) = 7.57 
Part of the uncertainty in this estimate comes from 
the range in possible Mbh values at nearly face-on 
inclinations (close to zero degrees), as shown in Fig- 
ure SJ Recent cross-correlation reverberation mapping 
resul ts quote statistical uncertainties of the order of 0.15 
dex jBentz et alJl2009b iDennev et aHl2010b iBarth et al.l 
120 llalTbh . but this neglects the uncertainty in the nor- 
malization factor, /, that is believe d to have an ob- 
ject to object scat ter of 0.44 dex (|Woo et al.l 120101: 
iGreene et al.l l2010bf ). Thus, our uncertainty in M B h 
for Mrk 50 is smaller that achieved by traditional rever- 
beration mapping estimates. Our independent measure- 
ment of Mbh can be used to estimate the appropriate 
value of f for Mrk 50 by com paring it to the virial esti- 
mate bv IBarth et~aTI (|2011bl) . M vir = / v 2 ct/G, where 
t and v are obtained from the cross-correlation of 
the continuum and broad line light curves and from 
the width of the broad line, respectively. We find 
l°gio / = 0.78lo'27i wn i cn agrees to within the errors 
with the commonly used mean values of log 10 (/) = 
0.74±g:i? from lOnken eTail (f200l . log 10 (/) = 0.72±g;?§ 
from IWoo etHI (|2010h . and log 10 (/) = 0.45±g g| from 
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iGraham et al.l (1201 ll ). We have used (/) to denote a nor- 
malization factor derived from large samples of reverber- 
ation mapped AGN A/bh estimates as distinct from the / 
value we measure individually for Mrk 50. A sample of 10 
independent black hole mass and / measurements with 
comparable uncertainties to Mrk 50 and Arp 151 would 
allow us to calculate a mean / value to ~ 0.3/vlO ~ 0.1 
dex uncertainty and to distinguish between these com- 
monly used mean values. 

An additional interesting feature of Figure|4]is the com- 
plex structure in the joint posterior distribution for the 
inclination angle and opening angle, a feature that was 
not seen in Arp 151. The joint posterior appears to have 
two distinct families of solutions, although one has al- 
most four times as much weight as the other. In an at- 
tempt to understand the origin of this structure, we sep- 
arated the posterior samples in the two modes in order 
to test whether they are correlated with any other pa- 
rameters (such as the inflow fraction), however we were 
unable to find any such correlations. Future improve- 
ments to the flexibility and realism of the model may 
enable us to rule out one of these modes, and hence con- 
strain the parameters more tightly and further reduce 
the uncertainties. 

While Mbh is well constrained, there are many ways to 
successfully model the large-scale structure of the rever- 
beration mapping data. This is illustrated by the degen- 
eracies in the posterior distributions plotted in Figure 
The quality of the model fits to the data are illustrated 
in Figure [21 including six model integrated H/3 flux light 
curves, an example of a model dataset of spectra for each 
epoch in the light curve, and two data spectra with the 
model spectra ovcrplotted. The smoothness of the mod- 
els compared to the data is illustrated in the spectral 
datasets of the data and model shown in the top panels of 
Figure [2] The Mrk 50 H/3 spectral profile did not change 
in shape drastically over the course of the LAMP 2011 
reverberation mapping campaign, and the model spec- 
tral profile is likewise very similar for all epochs. Even 
though the shape of the individual spectral profiles can 
be well-modeled, more sophisticated models will be re- 
quired to match the detail of the small-scale variability 
of the integrated H/3 data light curve. 

Note that the uncertainties quoted throughout this pa- 
per are determined from a Monte Carlo method, and are 
therefore subject to error themselves. As we are inter- 
ested in reducing the uncertainties on black hole mass es- 
timates from reverberation mapping data, it is important 
to quantify the uncertainty on the uncertainty. To inves- 
tigate this, we estimated the effective number of inde- 
pendent samples produced by our Diffusive Nested Sam- 
pling runs, by counting the number of times the explor- 
ing particles returned to the prior (allowing large steps 
to be taken) before climbing the likelihood peak again. 
Our effective number of independent samples was found 
to be ~ 180. We then generated samples of size 180 
from our full posterior sample, and determined the scat- 
ter in the resulting log 10 (A/BH) uncertainties to be 0.02. 
Thus, the uncertainty on the black hole mass for Mrk50 
is toir ± 0.02 dex. 

Previous attempts to understand the geometry and dy- 
namics of the BLR have focused on reconstructing the 
VRTF dKollatschnv fc Bischofl [200l iBentz etalll2010T: 



iDennev et al.|[2010f ). In the interests of comparing fu- 
ture transfer function studies to our physically motivated 
model of the BLR, we show three inferred VRTFs for Mrk 
50 in Figure [6] These three transfer functions were cho- 
sen out of the many inferred possible models for Mrk 50 
to show some of the variety in allowed transfer function 
shapes. The top left VRTF has a fairly typical shape 
and level of asymmetry, while the top right VRTF is 
more asymmetric than average. One measurement of 
the VRTF asymmetry is to compare the integral of the 
mean lag per velocity bin on either side of line center, 
corresponding to the zero velocity point in the middle 
right panel of Figure [6] By this measurement of asym- 
metry, 43% of the possible models inferred for Mrk 50 
have VRTFs that are less asymmetric than the top left 
VRTF, while only 8% of the possible models have VRTFs 
that are more asymmetric than the top right VRTF. 
The middle left transfer function illustrates the extent 
to which our inferred model for Mrk 50 can agree with 
the velocity - resolve d cross-correlation measurements by 
iBarth et al.l (|2011bD , shown by red errorbars in the mid- 
dle right panel of Figure [H This VRTF has the smallest 
Y 2 distance from th e cross-correlation measurements by 
IBarth et al.l (|201 lbf ) and models with this level of agree- 
ment (or better) have a probability of ~ 0.3%. The av- 
erage shape of the VRTF is also sh own in Figure El with 
the same velocity bins as used by IBarth et al.l ()2011bl ) 
for their cross-correlation based measurement. This av- 
erage VRTF is fairly symmetric, but the higher velocity 
bins have larger errorbars as a result of averaging over 
transfer functions that have asymmetries from either net 
inflowing or outflowing gas (see the dashed line in the 
middle right panel of Figure [6]) . 

Note that the average VRTF we infer and the re- 
sults obtain e d from cross-correlation measurements by 
IBarth et al.l (|2011b| ) do not all agree to within the l-a 
error bars. In order to understand the differences be- 
tween the time-lags as measured in our dynamical model 
and those measured through the cross-correlation proce- 
dure, we consider the ideal continuous noise-free case. In 
this case, the cross-correlation function (CCF) between 
the line and continuum light curves is the transfer func- 
tion convolved with the autocorrelation function (ACF) 
of the continuum light curve, which is the CCF of the 
light curve with itself. While the ACF is symmetric, 
the transfer function may be asymmetric, as we find for 
Mrk 50, so the CCF may also be asymmetric. One mea- 
surement of the cross-correlation time-lag often used to 
measure black hole mass is the CCF-weightcd mean lag, 
Teen j which is by definition affected by the asymmetry 
in the CCF. Therefore, in the case of asymmetric trans- 
fer functions, r ccn may not correspond to the mean lag 
of our dynamical model of the BLR. For the non-ideal 
direct comparison between cross-correlation mea- 
surements and the results of our dynamical modeling ap- 
proach is not straightforward, since the peak (or mean) 
of the CCF does not measure the true mean lag but only 
a noisy version of the convolution between the ACF and 
the transfer function. 

We explored this issue by run ning the cros s- correla tion 
technique as implemented by IBarth et al.l ()2011bl ) on 
light curves generated by models drawn from the pos- 
terior PDF for Mrk 50. For simplicity we considered 
noise-free light curves sampled in the same way as our 
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data. We find that the peak and CCF-wcightcd mean 
(^pcak and T ccn ) of the CCF can be systematically off by 
~ 1 — 2 light days with respect to the true mean lag of 
the model. The amount of the offset varies as a function 
of the actual shape of the transfer function as well as 
the details of the implementation of the cross-correlation 
algorithm. Thus, it is not surprising that we find system- 
atic differences of this order between our estimates of the 
mean lag and T ccn . Clarifying and quantifying systemat- 
ically the relationship between these two approaches as 
a function of BLR structure and data quality is an im- 
portant topic that goes beyond the scope of this paper 
and is left for future work. 

In conclusion, the analysis presented here provides new 
and unique insights into the geometry and kinematics of 
the BLR, and a Mbh estimate that is competitive with 
the most accurate methods. However, since our mod- 
eling uncertainties are greater than data uncertainties, 
more physical models that take into account the complex 
processes occurring in the BLR should allow for even bet- 
ter constraints. In the future, we plan to develop such 
models and apply them to large samples of reverberation 
mapping data. 
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